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Abstract 

In this article an hierarchy of terms in the expansion of the mul- 
tipolar Hamiltonian for non relativistic moving atoms in QED is con- 
sidered. The particular case of neutral composite systems of 2,3 and 4 
particles is considered. The proposed hierarchy is based on a scaling 
analysis of the multipolar Hamiltonian together with a multipolar ex- 
pansion of the coupling terms. We give explicit results up to the fourth 
order in the fine structure constant. 



1 Introduction 

The theory of multipolar Hamiltonians for atoms in QED has been mainly 
developped through the Power-Zienau-Wooley (PZW) transformation for 
the minimal couphng Hamiltonian. This formulation is well suited for the 
description of a confined system of charges and is particularly relevant for 
a neutral system where the vector potential disappears. For moving atoms 
there have been few proposals of such Hamiltonians [H, LBBL,GR] owing 
to the fact that, up to the first order in the fine structure a, symetries of 
the p ■ A and d ■ E schemes are the same. Higher order terms has been dis- 
cussed mostly in the framework of two level systems [LBBL, W], and more 



1 



recently for real atoms [BBB]. The analysis of usually disregarded terms 
[GR] appeals for a systematic proccdiirc that can lead to an unambiguous 
hierarchy of the terms in the multipolar expansions in order to perform a 
relevant perturbation theory. In this note we focus on QED of moving non 
relativistic light atoms, without external fields. These atoms are composite 
systems of two, three or four particles in order to make explicit the role of 
the increasing complexity of the problem. The point is to examine how the 
scaling technique performed from the beginning can be associated with mul- 
tipolar expansions in order to obtain a straightforward hierarchy of terms in 
the Hamiltonian. The second section reviews the QED multipolar Hamilto- 
nian and the relevant physical quantities, the third section presents a sketch 
of the scaling technique and proposes a way to determine which terms are 
predominant in the Hamiltonian and the fourth section makes a direct use 
of that choice within the multipolar expansion of the interaction terms and 
gives also the corresponding expansion for the minimal coupling Hamilto- 
nian. Finally, in the fifth section, the use of Jacobi coordinates is introduced 
and explicit expansions are given. 
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2 The multipolar Hamiltonian 

We consider a neutral system of iV charges {ea)i<a<N, with masses {ma)i<a<N, 
spins {Sa)i<a<N and Lande coefficient {ga)i<a<N- The PZW Hamiltonian 
for QED in Coulomb gauge is given by (see [GR]): 

iV - / ^ \ 2 



a=l " ^ a<b 



4:Tr£o\ra - n\ 



^ 1 ,(n(r) + P(r))2 5(r)2 



^Ma.Bira) + l f d\ (- 
a=l •' 



£^0 A*o 



+ ^^), (1) 



with 



p^ = -ifiS/r^, Ma = ga^r^Sa, C^EqUq = I (2) 



where 



A — 1)2 

-a^(A;)(AAe,(fc))e-^-) (3) 
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3 , (aA(A;)eA(fc)e^'^--al(fc)eA(A:)e-*'-) (4) 



and E{r) = Ii{r). (5) 

So 

Here a;(A;) = c\k\ and £A(fc)) = Ij 2, are real polarization vectors satisfying 

ex{k)-e^{k) =5xf,, k-£x{k) = 0, A,// = 1,2 (6) 

ax{k), a*^{k) are the usual annihilation and creation operators obeying the 
canonical commutation relations 

[ai{k),al{k)] = 0, [ax{k),a;{k')] = Sx^k - k') (7) 

where a" = a or a*. 

The polarization field relative to the center-of-mass coordinate R is: 

N 1 
a=l 

The magnetization fields relative to the center-of-mass coordinate R, for 
1 <a< N are 



N 1 

ea{r) = Y,eb{rb-R) / dA(A5„6 - ^(A - 1)) (^(r - i?i - A(r(, - (9) 
6=1 

where = E^i "i-a^'a with M = X^^i rria- 



3 Scaling 

3.1 General Remarks 

The scaling is a change of units which exhibits the pertubative character of 
the different terms in the Hamiltonian that we are concerned with. Here 
we use dilatation of particles coordinates and photon momenta indepen- 
dently. Then we have to establish which terms are to be predominant in the 
Hamiltonian of this composite (particles and photons) system. These terms 
will form the zero order Hamiltonian for the perturbative expansion. In a 
first step, we have to determine a relevant truncated zero order Hamiltonian 
extending to moving atoms the properties known for non moving atoms. 

In the hydrogcnoid model of atoms with independent electrons, typical 
binding energy for each electron is of the order of ^mc^(Za)^ [M] where 

2 

a = .Jt is the fine structure constant, Z the nucleus charge and m the 
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electron mass. This quantity, whicli neglects all the further corrections from 
inter electron repulsion sceening effects and exchange process, is an upper 
value for the total ionisation of an atom. Indeed the maximum atomic value 
for the first ionisation is that of the helium atom: 24.58 eV, of the order 
of (^mc^a^Z^)/2. As the multipolar formalism fits at best with confined 
system of charges we will disregard energies exceeding the first ionisation 
limit. Thus it is consistent to built a model where the Coulombic part, 
kinetic energy of the electrons and the free field energy have to be retained 
within the leading terms [FFG] . For the nuclear motion the question of the 
scaling factor arises: should we choose the same for nuclei and electrons? 
One point of view is to respect the consistency between different choices of 
coordinate systems for N particles systems: the natural one {ra}a=[i n] ^^'^ 
the Jacobi ones {Ra}a=[i n] which are generated by a definite building up of 
the center of mass and relative vectors and associated coordinates. As the 
constitutive set of equations corresponding to the transformation from one 
system to the other is a linear one, we have chosen to adopt a same scaling 
constant for nuclei and electrons in order to fit with this property. 

If the nucleus is the particle a = 1, the unperturbed zero order Hamil- 
tonian is the following one: 

where N = Z + 1. 
3.2 Dilatation 

We dilate the particles coordinates and photon momenta independently, 
i-e-, i{ra}a=[i,N],k) ({^}a=[i,N],r]ky M > 0, r? > 0, where /j, and r/ are 
dimensionless constants. This scaling induces a unitary transformation in 
the space of physical states that we now compute. 

The dilation k —>■ rjk induces an unitary operator in L'^{]R^) denoted by 

r, : 

(r^/)(fe) = 77^/2/(r/fc), r/>0, kelR" (11) 

where f e L'^{M^). 

We have for the inverse transformation 

{T-'9){k)=v-'/'9{^^) (12) 

Let J^ph be the Fock space for transversal photons. 

oo 

We have J'ph = L\1R^ ,(r'^)^" , where L'^{IR^,W'^)< = W. Here L'^ {M^ ,(r^)^^ 

n=0 
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is the symmetrized n-fold tensor product of {IB? ,W'^) . The dilation k 
rjk, T] > 0, induces an unitary operator on J^ph that we stih denote by F^. 
We have for tp = (V'^''^)n>o G 

(r^^)"(fci,Ai;-- - ;A;„,A„) = r?i"V^")(r?fci,Ai;-- - ;??A;„,A„), (13) 

kj G R^, Xj = 1,2, j = !,■■■ n. 

Thus we have 

r-'ax{k)Tri = v^/\x{vk) (14) 
r-'al{k)T, = r,'/'alivk) (15) 

because a^(fe) is the adjoint of ax{k). 

The right mathematical definition of the energy operator for the photons 
is given by 

1 f ,3 n(r)2 5(r)2 

- (Pr: — ^ H 5^ : (16) 

2 7 £o Mo 

where : : is the Wick normal order (sec [S]). 

Let Hf denote the energy operator for the photons. We then have : 

Hf = \Jd'r: + := ^ d'k \k\alik)axik). (17) 

A— 1,2 

Thus using (3), (4), (5), (19) and (14), (15) one gets: 



r,,B{r)T-' = V^B{rir), (18) 

r^n(r)r-i = r/2n(7?r), (19) 

Tr,E{r)r-^ =r]^E{7]r), (20) 

^vHfT-' = vHf. (21) 

Moreover 

r;V..r^ = /^v.„. (22) 

Furthermore, for any function V{ra), we have 

r;V(r„)r^ = y(^). (23) 



We have to transform the multipolar Hamiltonian by the total dilatation 
r acting both on particles coordinates and photon momenta: T = (F j_)^r^. 

Thus we now get: 
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rs(r)r-^ 
rn(r)r-i = r]^u{ 



•qr. 



THfT-' = vHf, 

3.3 Determination of the common scaling constant 

We first examine the action of F on the truncated Hamiltonian: 



(24a) 
(24b) 

(24c) 

(24d) 
(24e) 

(24f) 



FHoT-^ = 



N N 



.0=1 



a=2 



47r£oFa - ri 



n(r)2 B{r) 
£0 



a=l a=2 ' 



"(25) 
2 S(r)^ 



1^0 



Choosing the kinetic energy of the electrons, tlic free energy of the pho- 
tons and the Coulombic part as the leading terms of the same order we must 
have: 



/X = Zafi = r] 
We then get: 



V 



rj = (Za) , ji = Za and — = Za 



and we have: 



(26) 



(27) 



N N 



.a=l 



a=2 



he If-, Mr)"^ B(r)'^ 
r + - / d'^r 

ka-n 2 J 



£0 



+ 



(28) 



where the electronic ground state energy corresponding to the term within 
brackets is ^mc^Z as nia = m for a= [2, N]. 
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4 Scaled Multipolar expansion 



4.1 Scaled multipolar Hamiltonian 

The scaling conditions (24a — 24/, 27) and the equations (8 — 9) will now 
allow us to dilate all the terms of the multipolar Hamiltonian given by 
equation (1), once the expansion of squared bracket terms is performed. 



^ 1 2 ^ ^e' 1 Mr? B{r)\ 

JV TV 

+ Yl —P<^- / ^^'^ ®«(^) ^ ^(^) - E ^a-^('^a) (29) 

a=l a=l 

N 

+ E ^(/ 0«(^) ^ ^(^) + 2i^ / 



Constrains (26) and (27) fix the following behaviour for the involved 
operators: 



VHfV-^ = ^i^Hf, (30) 
VB{r)V-^ = n^B{ijLr), (31) 

rn(r)r-^ = /i^n(//r), (32) 

VE{r)V-'^ = i/E{iir). (33) 

Now it remains to compute the effect of the scaling transformation on 
terms including P(r) and 6a (r) : 
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AT 

.3 



V ea(r„ - i?). / dA / -n(r) 5{r - ii[R - X{ra - R)]) 
^ Jo J ^0 

.1 N 

= IJ? dXy^ eaiVa - ii).-n(/x[i? - \{Ta - i?)]), (34) 

r / dVe„(r)AS(r)r-i = /x3^e6(rfe-E)A / (A5„;, - ^(A - 1)) 
y 5(r) (^(r - i?i - \{r^ - R)) 



.1 TV 



= li^ I dX^{\6ab-^{X-l))ebin-R)A B{^i[R-\in-R)]). 

6=1 

(35) 

Finally one obtains for the multipolar Hamiltonian the following scaling: 



V (^^ + 



/t^^ „Jra-ri\ 2 J 

a=l a=2 ' ' 



3^ (IW ^ ^(0 ) 
N 

-/i2^M„.5(/xr„) (36) 

a=l 

a, 6=1 

(p..((ri,--R)AB(M|/;-A(rt-ii)|))) 

+ /Ei(f<'ABA..-^(A-l))e. 

a=l " •^'^ 6=1 

(rb-i?)A5(//[i2-A(rfe-i?)]))2 

As usual the vacuum polarization term, which scales as iJ,~^ has been for- 
getten in (36), and the last equality defines 
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4.2 Multipolar expansions 

The multipolar expansion of the expression (36) is straightforward. It relies 
on the usual Taylor series expansion for F = B or 11: 



where X stands for (r — Rg) or A(ra — R) . 

Note that \ra — R\ is of the order of |/s|ai(Z) where ai{Z) = mcZa/h. 
As the photons energy has to be lower to the value mc? which is used for 
the ultraviolet cutoff [BFS], it then results that \k\ < mc/H. Therefore the 
/i \ k.{ra — R) \ factor in the exponential terms of the field expansions (3) and 
(4) can be assumed to be small (mod 27r). 

One then gets the following explicit /x— expansions for these terms. 
Electric multipole terms: 




(37) 



n=0 





N 



(38) 



oo 



n=0 



Spin mutipole terms 



N 



a=l 




(39) 



oo 



(m) 



n=0 
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Magnetic multipole terms 

N 



-"J a,l3=l " 

{pa. {{rp -R)A B{i^[R - X{rp - R)]))) 

n=0 a,/3=l ^ 

oo 



n=0 

Diamagnetic multipole terms 

N , N 



a=l /3=1 

(r;3-i2)AS(/x[i?-A(r;3-i?)]))' 

A'' -I „i TV oo 



, 2mQ: /n ^ „ 
a=l •'^ /3=in=0 



M"A"(A5„^ - ^(A - l))e^ ((r^ - i?) A ^-^^^^^^Bir)^,^,^) }' 

(41) 

oo 

= E^"'"'^mm(/^) 

n=0 

This allows to perform direct expansion up to the required order in 
jj, = Za of the multipolar Hamiltonian for an atom considered as a confined 
neutral system of charges. Assuming that T'^,Tg ^T^^^Tj^j^ = for n < 
0, one gets for the interaction Hamiltonian the following multipolar 
expansion: 

^ oo 

-^THi^^T-^ = E {r^-\lA + rr'(/^) + Tlf\n) + r«-^(/x)} (42) 

^ n=0 

4.3 Multipolar expansion of the minimal coupling Hamilto- 
nian. 

For sake of completeness let us remind the expression of the minimal coupling 
Hamiltonian corresponding to the same system of charges as in (1): 
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^ 1 



^ 1 ^47r£oka -t"/?! 

a=l o</5 

M„.S(r.) + \jd\{^ + (43) 



where 



;^'tT;2"' \2eo(27r) a;(A;)/ 

(44) 



One easily computes that rA{r)T~^ = ii^A{iJ,r). For the other terms 
the preceding analysis remains inchanged. As a result one gets a scaled 
expression for the minimal coupling Hamiltonian: 



1 ^ ^ e 

+ 7 5Z I r ~ ^ 5Z —Pa-AifJ-ra) 

2<a</3' a=l 
AT AT 2 

- /x^ ^ M,.5(/x r„) + /x^ ^ 2^^(^^")' 

a=l a=l 

= Iri/or-^ + A<cr-\ (45) 

where the right hand side defines -ff^c- 

The multipolar expansion is the consequence of (37) for F = ^4 and 
X = (T-Ra): 



N 
a=l 

= E /^"^'(-l)'^ E ^;.,. ((^"- ^)-^-)" ^(.)|.^,^ (46) 
n=0 a=l " 

oo 

= EM"+'r^(M)- 

n=0 
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N 2 °° 
a=l " n=0 



Z^Z^2m« (n-l)\ ^^^''l-=M«- Z! ^l^Jk=Mi? ^^^^ 

i=0 a=l " ^ ^ 



00 



n=0 

Assuming that T^, T^^ = for n < 0, one gets for the minimal couphng 
interaction Hamiltonian the following expansion: 

-. 00 

-^ri/4,r-i = ^ {rr\^^) + Tr\i^) + ti-/{„)] (48) 

^ n=0 
Thus the correspondence between terms of the same order in ^ for H^p 
and Hmc can be established by direct comparison of the expressions (36) and 
(43) or (42) and (48) with multipolar expansions given by (38,39,40,41,46,47). 



5 Explicit Expansions for atoms with 1,2,3 elec- 
trons 

5.1 Jacobi vectors 

As the center of mass vector R and the {fa}i<:_a<N '^'''^ linearly dependent, 
one has to consider Jacobi vectors {Ra]i<a<N order to use independent 
variables. They are many such families of vectors, each one being related to a 
definite building of the center of mass which has to preserve the quantity / = 

Ea"^«''a = Y^aMaRl, whcrC (Ea "^«) ^1 = Y.a'^a^a^ ^1 = (Ea "^a) 

and Ma is the effective mass corresponding to the associated Ra- Usually 
one proceeds imaging what is done for two particles, pairing center of mass 
of (effective) particles getting at each step the center of mass, the relative 
vector and the associated masses. 



5.2 Octupole approximation for a neutral two body system 

We consider here the hydrogen atom. Here, ei = e, 62 = — e. We introduce 
the Jacobi vectors {Ra}a=i,2- 



We set 



111 „ rriiri + m2r2 

R2 = ri-r2, — = — + — ' -^1 = ^7 ■ 

M2 mi m2 M 

Ml = mi + m2, Pi = —itNu^, P2 = -iKVr^. 



12 



We obtain 

in = Ri + ^R2. [pi = P2 + ^Pi, 

In this new system of coordinates, H^p is given by: 

Hmp = Ho + We, (49) 

where 

p2 p2 2 
^ _ , _[2 ^ , 

° 2Mi 2M2 47r£o|i?2| ^' 



We = Hmp — Hq. 



1 



Thus, letting ^ = a in (38 — 42) the octupole approximation of —THmpT ^ 

is given by: Hq + W^, 
where 

and 

K = f'K + f^'^K + A^'^M + (51) 

with 

VF^i = -rg(/i) = ~eR2.E{aRi), 

w^2 = _rO(^) + ro,(^) + rl(^) 

= -(Ml + M2).B{aRi) + -^[Pi.eR2 A B(ai?i) + /i.c] 

1 4M9 1 

+ ^(1 - -j^)nP2.eR2 A i?(ai?i) + /i.c] 

1 4Mo 1 



w^' = ri^(/x) + r|(M) 



1 4M9 1 

(1 - lCr)'[^i-«^2 A (i?2.V.)i?(r)|_^^ + /i.c] 



4Mr Ml ' ' . V . ,/ V /|.=„fli 
+ 6^2 " ^^^''^•''^^ (^2.V.)S(r)|_,^ + ^.c] 
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< = r|^(/i) + T|(/x) + TOfAf(/x) 

= ^— (l-^)[Pi.ei?2 A(i?2.V,)^5(r)| , +h.c] 

12Ml Ml ■'|r = <.Hi J 

Let us remark that, at the second order in fi, the term i?2.(-R2-Vr)-E(r)|^^^^^ 
does not appear in [GR] because the authors only consider the dipolar ap- 
proximation. The quadrupolar electric term appears as soon as one considers 
the quadrupolar approximation. 



5.3 Quadrupole-approximation for a neutral three-body sys- 
tem 

Now we consider a neutral 3 body system as, for example, the helium atom 
He'^. Let be the Hamiltonian of the considered system. 

Consider a system formed of 3 particles Ai,A2, A3 of respective masses 
mi,m2, ms, with charges ei, 62, 63, of positions ri, r2, rs and of impulsion 

Pi, P2, P3- 

We introduce the Jacobi vectors associated to the following partition: 
a = {(l,2), 3}. Itis the case of a neutral atom formed of a nucleus of charge 
— 2e and of two electrons each one having a charge e. 

We set 



mm + m2r2 + msrs 
R2 = r2 — ri, iti = — , Ml = 1711+1712 + ms 



R3 = 
Pi = 
We obtain 



^3 



miri m2r2 



M2 

ihVR,, P2 



m 



Ml 

M2 = mi + m2- 



-ihV 



R2, 



-ihV 



( 

ri 


= Ri- 


Ml ^3 






Pi = 


Mi-^1 


-^^3-^2 


i r2 


= Ri- 


Mr^3 + 






P2 = 




-^^3+P2 


J3 


= Ri + 


M2 p 






,P3 = 




+ P3 + 



In the case of the Helium atom we have mi = m2, &i 
63 = -2e. 

In this new system of coordinates, iJ^^ is given by: 



62 = e and 
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+ 



47r£o|i?2| ^TTEolRs + ^i?2| 47r£o|-R3 - 5P2I 

3 

Ma.B{ra) + Hf - / E{r).P{r)d^r. 
^1 -^^^ 



We write iJ^^ under the following form: 



with 



^He Pi Pi Pi 2e' 



2Mi 2^12,3 2mi 47r£o|ii:2| 47reo|-R3 + 5^2! 

2e^ ^ 
47r£o|r3-^i?2| ^ 



and 



Here ^- = ^5^- + 

/ii2,3 2mi ms 



Letting A?^ = 3 in (38 — 42) we get at the second order in /x = 2a: 
where, 

° " 2Mi + 2^i2,3^mi + 2|i?2| \R3 + lR2\ \R3-lR2\'^ ^' 
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with 



K = -'^e(i^) = 2ei?3.^(2ai?i), 
W^^=-Tl{fi) + TUl^)+Tk{fi) 

= -(Ml + M2 + M3).B{2aRi) - -^[Pi.eRs A B{2aRi) + h.c] 

Ml 

+ - %^)^ [P3-ei?3 A B(2ai?i) + h.c] 



2^12,3 
1 



Ml 



(58) 



+ [P2.eR2 A B{2aRi) + /i.c] 

4mi 

- e[(l - ^)^i?3.(i?3.V.) + ii?2.(i22.V.)]£;(r)i^^,^,^ . 

5.4 Quadrupole-approximation for a neutral four-body sys- 
tem 

We are finally interested in the case of neutral moving system formed of four 
bodies. For this system we consider the following partition for {1,2,3,4} : 

c = {a, 4}, a = {1,2, 3}. 

This partition correspond to the lithium atom. 
We define the Jacobi vectors as follows: 



R2 

Ri 
P2 



min + m2r2 lu u \ 

n - r2, -K3 = r3 — ;— , Ri = (-K4, R2,) 



mi + m2 



miTi + m2r2 + mars 



mi + m2 + ms 
miri + m2r2 + msrs + m4r4 



M 



-ihV 



R2' 



M4 = mi + 777-2 + 1^13 
, Ml = 7771 + ^^2 + 7773 + 7774 

p^ = -invR^, Pi = -ihVR,. 



We suppose here that 777 1 = 7772 

7774 = 7772. 

We then obtain 



7 , r-4 = —ifov R^, ri 

7773 and we set 7771 = 7772 = 7773 and 





= Ri- 


Ml -"-4 


|i?3 + 


r2 


= Ri- 


Ml -"4 


\Rz - 




= Ri- 


Ml -"4 


iRz 




= Ri + 


Ml -"4 





\R2 



Pi = 


Mi^J^ 


iP4 - ^^3 + P2 


^ P2 = 


Mi-^l 


|P4 - 5P3 - P2 


P3 = 


Ml 1 


|P4 + P3 


P4 = 


Mi-^l ^ 


P4. 
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We recalls that the Hamiltonian i^^p of the considered system is given by: 



mp 

4 



a=l " l<o</3<4 ^' 



Q;=l 

We now consider that ei = 62 = £3 = e and 64 = — 3e, the lithium atom. 

In the new system of coordinates the Hamiltonian -ff^p is written as 
follows: 

■(^Pi+J'4+ / <i='re4(r)AB(r))' 



1 ^mi 
2mi 

2m2 



+ ^ r^ + - + 



ATTeo\R2\ 47r£o|-R3 - 5^1:21 47r£o|i23 + 
3e2 3e2 



47r£o|i?4 + ^Rs - \R2\ 47r£o|-R4 + \Rz + 5-R2I 

+ Hf-Y, Mo^.Bira) - j E{r).P{r)d\. 



3e2 



47r£o|i?4 - l^sl 
We then write under the following form: 



where 



2M\ 6mim2 4mi mi 

+ ^ 77^ + -. + 



47r£o|i22| 47r£o|it;3 - 5i?2| 47r£o|ii:3 + ^ii;2| 

3e2 3e2 



(59) 



47r£o|i?4 + I-R3 - 5^2! 47r£o|i?4 + 5-^3 + \R2\ 
3e2 



47r£o|i?4 - 5-^3 1 



2^ +^/' 
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and 



We = ~Yl ^oc-B{Ri) - I E{r).P{r)d\ 

4 4 

+ 2^(j'i-/ d\{Y.e^{T)^B{r) + J dV(^e,(r)AS(r).Pi) 

3=1 j=l 
3 3 

+ ^(P4.y" o?Ve4(r) AB(r) + y dVe4(r) AB(r).P4) 
-^(Pa.y" rfV(ei(r) + 62(0-263(0) A 5(r) (60) 
+ y" dV(ei(r) + 62(0 - 2e3(r)) A ^(O.Pa) 
+ ^(P2.y" c?V(6i(r)-62(r)) AP(r) 
+ J d^riei{r)-Q2ir))AB{r).P2) 

+ 2^{( y dVei(r)AP(r))' + ( ^ c/Ve2(r) A B(r))'} 

+ -^( / dV63(r) A5(r))^ + -^( / dV64(r) AS(r))^ 
zmi J Zm2 J 



Let 

= ^^Hi^v^~'- (61) 
Setting A/" = 4 in (38 — 42) we obtain at the second order in /x = 3a: 

Hj:' = H^' + iiWl + i?Wl (62) 

with 



Li _ Pi I -^1 p 2 I 3 2 , 1 p 2 
— TTirr + ^ + -j -rs H -^2 

/Ml bmim2 4mi mi 

/ic he 

^ W+W^\ ^^^^ 

3hc She 3hc 

I-R4 + \Rz — \R2\ \Ra + \Ri + ^-^2! IP4 — 
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Wl = -T^iij) = -3eR4.E{3aRi), 

= -[Ml + Ms + Ms + M4].B{3aRi) 

^('^)|.=3a«i 

+ -^[Pi.d A B{3aRi) + h.c] + —[P^.d A B{3aRi) + h.c^ 
1M\ 12 

1 1 

+ [Ps.ei^s A Bi3aRi) + h.c] + \P2.eR2 A B(3aRi) + h.c]. 

4mi 4mi 

(64) 

We observe that the several approximations for the neutral composite 
systems of 2, 3 and 4 particles we have obtained are exactly the same as the 
ones got by [F]. In [F] the author first computes the multipolar expansion 
and then implements the scaling analysis. 



6 Conclusion: 

We have found again (sec [GR]) the electric dipolar term at first order in 
jj, and at second order in fj, the magnetic dipolar term together with the 
Rontgen one that already appear in the dipolar approximation of Hmp- But 
now, considering the full multipolar expansion in /x, we have shown that 
an extra term, i.e., an electric quadrupolar term appears also at the second 
order in /j,. 

This strongly suggests that we must finally consider in energy computa- 
tion the following Hamiltonian instead of the usual dipolar one: 

H(^l = Ho + {-T% + {-T^ + Tl, + Th}). (65) 

As the symmetries of the involved terms are different, a new level of com- 
plexity has to be taken into account. For instance for the computation of 
the spontaneous emission of these atomic system, the deductions of [WB] 
have to be consider again. 
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